Primordial magnetic field limits from cosmological data 
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We study limits on a primordial magnetic field arising from cosmological data, including that 
from big bang nucleosynthesis, cosmic microwave background polarization plane Faraday rotation 
limits, and large-scale structure formation. We show that the physically-relevant quantity is the 
value of the effective magnetic field, and limits on it are independent of how the magnetic field was 
generated. 

PACS numbers: 98.70.Vc, 98.80.-k 
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I. INTRODUCTION 



There is much interest in the origin of the coherent 
part of the large-scale /iG magnetic fields in galaxies p]]. 1 
A leading possible explanation is that these large-scale 
magnetic fields are the amplified remnants of a primor- 
dial seed magnetic field generated in the early Universe 
[1-0| • Such early magnetogenesis models include those in 
which seed magnetic field generation occurs during infla- 
tion or shortly thereafter, or during a cosmological phase 
transition (such as the electroweak or QCD transition). 
Clearly the strength of the seed magnetic field should be 
small enough so as to not generate a larger than observed 
cosmological anisotropy. Magnetic field energy density 
contributes to the relativistic (radiation) energy density 
and thus another requirement is that it not exceed the 
big bang nucleosynthesis (BBN) bound on the radiation 
energy density. 

There are two main questions that need to be an- 
swered: (1) Are the amplitude and statistical proper- 
ties of any of these seed magnetic fields such that, af- 
ter amplification by a realistic model, they can explain 
the strengths and correlation lengths of the observed 
magnetic fields in large-scale structures (LSSs) such as 
galaxies? and, (2) Are any of the seed magnetic fields 
detectable through cosmological observations, such as 
cosmic microwave background (CMB) measurements 2 or 



LSS observations? And, if yes, what are the observational 
constraints on such primordial magnetic fields? 

In this paper we focus on the second question and 
consider two cosmological consequences of a primordial 
magnetic field, Faraday rotation of the CMB polariza- 
tion plane and the effect on LSS formation. As two of 
the effects of a primordial magnetic field, these have been 
widely discussed in the literature. See Refs. [H-[l3[ for dis- 
cussions of magnetic field induced CMB anisotropics, 3 
Refs. 15l-[r7| for the Faraday rotation effect, and Refs. 
[18h22| for effects of a primordial magnetic field on LSS 
formation. 

It has become conventional to derive the cosmological 
effects of a seed magnetic field by using a magnetic field 
spectral shape (parametrized by the spectral index ns) 
and the smoothed value of the magnetic field (B\) at a 
given scale A (which is usually taken to be 1 Mpc). We 
develop here a different and more correct formalism based 
on the effective magnetic field value that is determined 
by the total energy density of the magnetic field. As 
a striking consequence, we show that even an extremely 
small smoothed magnetic field of 1CU 29 G at 1 Mpc, with 
the Batchelor spectral shape (ng = 2) at large scales, can 
leave detectable signatures in CMB or LSS statistics. 4 

We also show that the conventional approach based 
on the smoothed magnetic field results in some confu- 
sion when considering phase-transition generated mag- 
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1 On larger scales there have been two recent claims of an observed 
lower limit of order 10 —15 G on the intergalactic magnetic field 
3, , as well as one claimed detection of a field of this strength 
4|. Prior to these observations, the intergalactic magnetic field 
was observationally only limited to be smaller than a few nG. 

2 A cosmological magnetic field induces all three kinds of gravita- 



tional perturbations, scalar, vector, and tensor; all three of which 
contribute to CMB temperature and polarization anisotropics. 

3 The effects of a homogeneous magnetic field on the scalar mode 
of CMB fluctuations, including the resulting non-Gaussianity of 
the CMB temperature map, are discussed in Refs. [l4|l , 

4 This strong limit on the primordial magnetic field is the conse- 
quence of the BBN bound and the sharp shape of the magnetic 
field at large scales [23l. |24| . The low efficiency of gravitational 
wave production by a cosmological magnetic field |24| results in 
a weaker limit on the seed field from the direct detection of the 
induced gravitational wave signal [25t . 
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netic fields [26[ with spectral shape sharper (on large 
scales) than the white noise spectrum (i.e. with ub > 0). 
In this case the total energy density of the magnetic field, 
Pb, is mainly concentrated on large wavenumbers (small 
length scales). In what follows we show that primor- 
dial magnetic field effects on cosmological scales are de- 
termined not by the amplitude of the magnetic field on 
these scales, but rather by the total energy density of the 
magnetic field. 

Here we consider limits on a seed magnetic field from 
the observational constraint on the CMB polarization 
plane rotation angle and observational constraints on the 
formation of the first bound structures. Both tests give 
comparable limits on the effective magnetic field, ranging 
from 10~ 9 G to 10~ 7 G, depending on the spectral shape 
of the magnetic field. Note that these limits are of order 
of the BBN bound. The best limit on the seed magnetic 
field is for the scale-invariant case that can be generated 
during inflation [27ll28||, 

In a Universe with only scalar mode perturbations the 
CMB B-polarization signal vanishes. The CMB B po- 
larization signal can arise from vector or tensor pertur- 
bations, and thus B-polarization detection based tests 
are powerful tools for probing non-standard cosmologi- 
cal models and the relic gravitational wave background. 
Since a cosmological magnetic field can source a CMB B- 
polarization signal, a crucial test to limit the magnetic 
field is based on CMB £?-polarization measurements. In 
a separate paper we will address the cosmological mag- 
netic field energy density limits that can result from this 
test. 

The structure of our paper is as follows. In Sec. II 
we review magnetic field statistical properties. In Sec. 
Ill we examine the CMB polarization plane Faraday ro- 
tation effect and the resulting magnetic field limits. In 
Sec. IV we consider the influence of a magnetic field on 
LSS statistics and determine the resulting limits on the 
magnetic field. We discuss our results and conclude in 
Sec. V. 



II. MODELING THE PRIMORDIAL 
MAGNETIC FIELD 

A stochastic Gaussian magnetic field is fully described 
by its two-point correlation function. For simplicity we 
consider here the case of the non-helical magnetic field for 
which the two-point correlation function in wavenumber 
space is [ll[ 



is the power spectrum of the magnetic field. 5 

We define the smoothed magnetic field B\ through the 
mean-square magnetic field [ill ] 



B, 



(B(x)-B(x))|; 



(3) 



where the smoothing is done on a comoving length A with 
a Gaussian smoothing kernel function oc exp[— x 2 /X 2 ]. 
Corresponding to the smoothing length A is the smooth- 
ing wavenumber k\ = 2tt/X. The power spectrum Pe(fc) 
is assumed to be a simple power law on large scales, 
k < kjj (where ko is the cutoff wavenumber), 



P B (k)=P BO k nE 



2tt 2 \ 3 B 2 
r(n B /2 + 3/2) 



(\ky 



(4) 



and assumed to vanish on small scales where k > hp. 
The energy density of the magnetic field is (l3j 



Pb(vo) 



B 2 (k D \) n °+ 3 
87rr(n s /2 + 5/2)' 



(5) 



We define the effective magnetic field B e g through ps = 
B 2 s /(8ir) and thus we get for the scale-invariant spec- 
trum (riB — —3 |27j) B e g = B\ for all values of A. The 
scale-invariant case is the only case where the values of 
the effective and smoothed fields coincide. 

We need to define the magnetic field cut-off wavenum- 
ber ko- We assume that the cut-off scale is determined 
by the Alfven wave damping scale kjy ~ vaLs where va 
is the Alfven velocity and Ls the Silk damping scale [l(| ■ 
Such a description is more appropriate when we are deal- 
ing with an homogeneous magnetic field and the Alfven 
waves are the fluctuations Bi(x) with respect to a back- 
ground homogeneous magnetic field B (|Bi| -C |B |). 
In the case of the stochastic magnetic field we generalize 
the Alfven velocity definition, see Ref. [tj| . by refer- 
ring to the analogy between the effective magnetic field 
and the homogeneous magnetic field. Assuming that the 
Alfven velocity is determined by B c g , a simple computa- 
tion gives the expression of ko in terms of B e g [24| : 



lMpc" 



= 1.4 4 



(2n) n ^+ 3 h /10" 7 G 



r(n B /2 + 5/2) V B cB 



(6) 



Here h is the Hubble constant in units of 100 km s" 1 
Mpc -1 . The BBN limit on the effective magnetic field 
strength, B cS < 8.4 x 10~ 7 G [Hjj, gives an upper limit 
on the cut-off wavenumber kp, 



( 9 \(n B +3)/2 

fcg BN > o.nh 1 ' 2 , — M PC - 

D - rV2( ns /2 + 5/2) 



(7) 



(5*(k)5,(k')) = (27r)W 3 >(k - k')P y (k)P B (fc). (1) 

Here i and j are spatial indices, i, j £ (1,2,3), ki — ki/k a 
unit wavevector, Py(k) = <5,j — kikj the transverse plane 
projector, <j( 3 )(k— k') the Dirac delta function, and Pe(fc) 



5 We use 

B 3 -(k) = I Ae^BjW, B J (x) = J e-^BjpO, (2) 

when Fourier transforming between position and wavenumber 
spaces. We assume flat spatial hypersurfaces (consistent with 
current observational indications, 129(1 '). 
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In the case of an extremely large magnetic field it is pos- 
sible to have A^> > 1 Mpc. At this point it would seem 
unreasonable (unjustified) to consider a smoothing scale 
A = 1 Mpc as is conventionally done. 



III. CMB POLARIZATION PLANE ROTATION 

The presence of a primordial magnetic field during re- 
combination causes a rotation of the CMB polarization 
plane through the Faraday effect [l5|. The rms rotation 
angle a rms = ((a 2 )) 1 / 2 induced by a stochastic magnetic 
field with smoothed amplitude B\ and spectral index ub 
is given by 



or 



E 



2/ + 1 

4tT 



CP 



(8) 



where the rotation multipole power spectrum CP is [Tj 



CP- 
91(1 + 1) 



(4^)3 g 2 ^r> B /2 + 3/2) 



n B +3 pxs 



(9) 

dxx nB jf(x). 



Here rjo is the conformal time today, vq is the CMB pho- 
ton frequency, q 2 = 1/137 is the squared elementary 
charge in cgs units, ji(x) is a Bessel function with ar- 
gument x — kr/Q, and xs — kgrjQ where ks — 2 Mpc -1 is 
the Silk damping scale. In the case of an extreme mag- 
netic field which just satisfies the BBN bound, kn might 
become less than the Silk damping scale. In this case 
the upper limit in the integral above must be replaced 
by x D = k D r/ . 

In terms of -B e ff, Eq. ® can be rewritten in the fol- 
lowing form, 



CP ~ 1.6 x 10 



1(1 + 1) 



x- 



n B 



(k D Vo) nB+3 
3 r xs 



B, 



cir 



incy v 



2 / 100 GHz x 4 



dxx nB jf(x), 



(10) 



and, as a result, 
a rmB ~ 0.14° 



B, 



off 



InG 



100 GHz 



x 



oo 

E 

.1=0 



[21 + 1)1(1 + 1] 



(fcr.r?o) ( " B+3)/2 

1/2 



dxx nB jf(x) 



.(11) 



It is of interest to compare Eq. (ITT1) with the correspond- 
ing result, Eq. (2) of Ref. [lj|, derived for a homogeneous 
magnetic field and at frequency vq — 30 GHz, 



1.6° 



InG 



^30 GH Z y 



(12) 



Both expressions agree for ng — > — 3 after accounting for 
J2i(2l + l)ji( x ) = 1 an d the fact that Bessel functions 
peak at x ~ I for given I (see Appendix A) . 




FIG. 1: Rms rotation angle function of spectral 

index n_e for the case when B c ff — 1 nG and vo = 100 GHz. 
Circles correspond to the computed values. 
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FIG. 2: Effective magnetic field limits set by the measure- 
ment of the rotation angle a rms for different spectral indices 
(n_B = —3,-2,-1,0,1,2, from bottom to top). The hori- 
zontal solid line shows the upper limit set by BBN. Vertical 
dashed lines correspond to the angles a rms = 3.16° that is set 
by the BBN limit on the effective magnetic field with spec- 
tral index n B = 2 and a rms = 4.4° set by the WMAP 7-year 
data. The numerical values of the effective magnetic field con- 
straints (in nG at 100 GHz) from the a rms = 4.4° limit are 
shown on the graph for each spectral index value. 



Figure [T] shows the rms rotation angle a rms , Eq. (jlll) . 
as a function of the spectral index ns when the effective 
magnetic field is normalized to be 10~ 9 G. The WMAP 
7-year data limits the rms rotation angle to be less then 
4.4° at 95% C L. [30|. This allows us to limit the effective 
magnetic field as shown in Fig. [5J 
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IV. LARGE-SCALE STRUCTURE 

A primordial tangled magnetic field can also induce the 
formation of structures in the Universe. In particular, 
these fields can play an important role in the formation 
of first structures (see, e.g. Refs. @, HHH H [H| ) . 

The magnetic-field-induced matter power spectrum 
P(k) is oc k 4 for n B > —1.5 and oc k 2nB+T for n B < —1.5 
[19l |32| . The cut-off scale of the power spectrum is deter- 
mined by the larger of the magnetic Jeans' wavenumbcr 
kj and the thermal Jeans' wavenumber fctherm (for a de- 
tailed discussion, see, e.g. Ref. [HI). Here the magnetic 
Jeans' wavenumber is (see, e.g. Ref. [32| ) 

kj ~ (230(" B+3 >/ 2 x 13.8) 2 /(" B + 5 ) ( ^] Mpc- 1 . 

\ B s J 

(13) 

Unlike the ACDM matter power spectrum, the magnetic- 
field-induced matter power spectrum increases at small 
scales and can exceed the ACDM matter one at small 
scales (for a comparison of these two spectra, see, e.g. 
Fig. 3 of Ref. [19]). And, therefore, one of the more 
important contributions of the additional power induced 
by magnetic fields is to the formation of the first struc- 
tures in the Universe (e.g. Refs. [18|, 12CJ, |2l| and references 
therein). 

In Fig.|3]we show the linear mass dispersion a(M) for 
matter power spectra induced by a primordial magnetic 
field with B c g = 6 nG at z = 10 for different values of Ub- 
Notable features of Fig.[3]are: (a) the mass dispersion on 
small scales is larger for a larger value of tib] and, (b) 
for riB > —1.5, the mass dispersion drops more sharply 
at larger scales than for ub < —1.5. We focus here on 
the mass dispersion on the smallest scales, as these scales 
are more relevant for the formation of the first structures 
in the Universe. These first structures were responsible 
for the reionization of the Universe at z ~ 10. To obtain 
meaningful constraints on i? c ff from the formation of first 
structures, we need to know how the curves shown in Fig. 
[3] vary as B e s is changed and as the Universe evolves. 

The mass dispersion a(M, z) evolves with the time 
dependence of the growing mode of the linear density 
perturbations sourced by the primordial magnetic field 
HHH^. The growing mode is oc a(t), the scale factor, at 
high redshifts, the same as in the "standard" ACDM case 
without a magnetic field. To account for this evolution 
the curves corresponding to a in Fig. [3] must be scaled 
by roughly a factor of ~ 11/(1 + z) for redshifts z 1. 

It can be shown that the value of a at the small- 
est scales (M ~ 1O 6 M ) is invariant under a change 
in B e g if the cut-off scale is determined by fcj: an in- 
crease/decrease in the value of i? e ff is compensated by a 
decrease/increase in the value of kj. However, if B c g is 
decreased to a value at which fcthcrm — ^ J > t nen the value 
of a decreases with a decrease in B g , as the cut-off scale 
becomes independent of the value of B c g. 

It has been shown that the dissipation of magnetic 
fields in the post-recombination era can substantially al- 




Mass (Solar Mass) 



FIG. 3: The mass dispersion at z = 10 for £> c ff = 6 nG as 
a function of magnetic field power spectral index tib- From 
top to bottom (at the left hand side of the plot), the curves 
correspond to nj = 2, 1, 0, —1, —2, —2.8. 
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FIG. 4: Constraint on the magnetic field strength B e ff as a 
function of the power spectral index ub- 



ter the thermal and ionization history of the universe 
[l~8l I20I . |22| . In particular, this dissipation raises the 
matter temperature and therefore the Jeans' scale in the 
IGM. For B c ff > 1 nG the matter temperature rises to 
~ 10 4 K as early as z > 100, 2Q], resulting in a steep rise 
in the Jeans' scale as compared to the usual case. The 
Jeans' wave number corresponding to this temperature 
is fcthorm ^ 10 Mpc" 1 (see, e.g. Fig. 4 of Ref. [22j|). 

WMAP results show that the Universe reionized at 
z ~ 10. This reionization was caused by the non-linear 
collapse of the first structures, followed by star forma- 
tion and the emission of UV photons from the collapsed 
halos. For a virialized structure in the spherical collapse 
model, the linear mass dispersion a ~ 1.7. This implies 
that the value of a at the scales of interest at z ~ 10 
is not expected to be much higher than 1.7. Consider 
the ub = 2 model in Fig. [3] the value of mass disper- 
sion at the smallest scales is ~ 100, which means that 
the first structures formed at z ~ 650 in this case (the 
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redshift of the collapse of first structures is ~ 6.5<r max , 
where cr max is the maximum value of a at z ~ 10), which 
can certainly be ruled out by the WMAP data on CMB 
anisotropics. A similar arguments can be used to rule out 
almost all the models shown in Fig. [3] Only the nearly 
scale-invariant models with tib — —3 do not put strong 
constraints on the strength of the magnetic field. As ar- 
gued above, the value of mass dispersion at the smallest 
scales to collapse is nearly independent of the magnetic 
field strength unless B c g decreases to a value such that 
kj = ^thcrm- In this case, the value of a decreases below 
those shown in Fig. [3J We have explored a wide range 
of -B e ff for the range of spectral indices shown in Fig. [3] 
We find that the range of acceptable values is 1-3 nG. In 
Fig. |4] we show the B e g corresponding to /cj = fc t herm- 
Notwithstanding various complications discussed above, 
this figure gives a rough sense of the acceptable range of 
-B e fT over the entire range of ub ■ 

In the foregoing, we neglect the impact of the ACDM 
model on the process of reionization. As the density 
fields induced by the ACDM model and the magnetic 
field are uncorrelated, the matter power spectra owing 
to these two physical phenomena would add in quadra- 
ture. The smallest structures to collapse at z ~ 10 in the 
WMAP-normalized ACDM model are 2.5a fluctuations 
of the density field as opposed to the magnetic field case 
where ler collapse is possible (Fig. 3). This means the 
number of collapsed halos is more abundant in the latter 
case. Therefore, depending on the star-formation history, 
if the magnctic-ficld-induced halo collapse made an im- 
portant contribution to the reionization process, the far 
rarer halos from ACDM would have made a negligible 
impact (for further details and references see Ref. [2l|). 



V. CONCLUSIONS 
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FIG. 5: The squared spherical Bessel functions jf(x) for I = 
5000 (top) and I = 10000 (bottom). Clearly jf(x) peaks at 
x ~ I. 
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FIG. 6: The sum of the squared spherical Bessel function 
Yli=o(?l + l)j?( x ) f° r lc ~ xs ^ 16000. The sum converges 
to 1 (horizontal solid line). 

bounds on the intergalactic magnetic field if the field 
was generated in the early Universe with spectral shape 
Ub < 1. This favors the inflationary magnetogenesis sce- 
nario. 



In this paper we study the large-scale imprints of a cos- 
mological magnetic field, such as the rotation of the CMB 
polarization plane and formation of the first bound struc- 
tures. We derive the corresponding limits on a primordial 
magnetic field energy density, expressed as limits on the 
effective value of the magnetic field, B e g. These limits 
are identical to limits on the smoothed magnetic field B\ 
(smoothed over a length scale A that is conventionally 
taken to be 1 Mpc) only in the case of the scale-invariant 
magnetic field (when ns = —3). For a steep magnetic 
field with spectral index ns — 2 the difference between 
B\=i Mpc and B e g is enormous (greater than 10 15 ). We 
show that using the smoothed magnetic field can result in 
some confusion; e.g. an extremely small smoothed mag- 
netic field on large scales does not mean that this field 
cannot leave observable traces on cosmological scales. 

An intergalactic magnetic field of effective value larger 
than 1-10 nG (with, depending on magnetic spectral 
index, corresponding values of B\=i Mpc in the range 
10~ 8 — 10~ 26 G) is ruled out by cosmological data. These 
limits of 1-10 nG are consistent with recent observational 
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Appendix A: Evaluating the right hand side of Eq. 
(fTTj) when n B -> -3 

The \JriB + 3 factor in the numerator of the right 
hand side of Eq. (fTTj) is compensated by a correspond- 
ing \j \/tib + 3 from the Bessel function integral when 
the spectral index ns ~ > — 3 and so the expression for 
a rms remains finite in this limit. To establish this we use 
properties of the Bessel function. Recall that jf(x) peaks 
at x ~ I, as shown in Fig. [5] This allows us to replace 
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the factor 1(1 + l)jf(x) by x 2 jf(x) (the accuracy of this 
approximation is of order 15-20%). The next step is to 
perform the sum over /. It is obvious that there is cut- 
off multipole number Iq that corresponds to the cut-off 
wavenumber, lc ~ min(xD, %s)- Now jf(x) satisfies 

oo 

£(2Z + I)j?(a0=l, (Al) 

1=0 



while we are interested in computing 5Zj= (2Z + l)jf{ x )- 
The Silk damping scale cutoff multipole number is Is — 
16000, 117 . Figure [6] shows that the sum to Is converges 
to 1. 
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